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We develop the two-electron attosecond streak camera under realistic conditions using a quasi- 
classical model. We assume extreme ultra-violet (XUV) attosecond pulses with a full width at half 
maximum (FWHM) of 24 attoseconds, centered at 120 eV and a streaking infrared laser field of 
1600 nm, and intensity of 1.8 x 10^^ W/cm^. The proposed method is shown to be capable to time 
resolve intra-atomic collisions in double ionization. 
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Time-resolving correlated electron processes is one of 
the driving forces behind the large scale effort to push the 
frontiers of attosecond science. Attosecond science offers 
time resolution through XUV-pulses. However, pump- 
probe experiments using attosecond pulses are techni- 
cally very challenging. Hence, the streaking of photo- 
electrons with an infrared (IR) laser field has become 
a successful technique for bringing time resolution to 
photo-ionization. The paradigmatic attosecond streak 
camera [il, 0, and its development into FROG CRAB 
(Frequency Resolved Optical Gating for Complete Re- 
construction of Attosecond Bursts) [3| originally aiming 
to characterize attosecond extreme ultra-violet (XUV)- 
pulses, have been the underlying concept for studies re- 
solving delayed time emission from atoms [1-[8| and solids 

i- 

To study the electron correlation in single photon dou- 
ble ionization we recently formulated the two-electron 
streak camera [13] • Specifically, it was shown that the 
intra-atomic knock-out process [llj can be associated 
with a delay, i.e. a time between photo-absorption and 
ejection of two electrons in the continuum. The delay is 
encoded in the inter-electronic angle of escape as a func- 
tion of the phase of the IR laser field. In addition, the 
two-electron streak camera can time-resolve delays cor- 
responding to different energy sharings, between the two 
electrons, and to different ionization mechanisms [l2j . 
However, previous work considered only discrete photon 
energies and instantaneous photo absorption [lo|, [HI • 

Here, we remove this severe limitation and extend the 
two-electron streak camera to realistic attosecond pulses. 
By resolving the bandwidth of an XUV-pulse in the sum- 
energy of two emitted electrons we construct the two- 
electron equivalent of FROG to obtain a complete picture 
of the single photon double ionization process. Specifi- 
cally, in FROG [l3| one extracts from a two-dimensional 
data set (FROG-trace) the complete characteristics of an 
optical pulse. In a similar manner in FROG- CRAB [sl 
one retrieves the spectral phases and amplitudes of an 



attosecond pulse. Here, we assume a fully characterized, 
transform limited XUV attosecond pulse, to obtain in- 
formation about the delay of two-electron emission after 
absorption of a photon from the attosecond pulse. We 
use as FROG-like trace for double ionization the inter- 
electronic angle of escape as a function of the phase of 
the IR laser field. A simple algorithm is devised for ex- 
tracting the intra-atomic two-electron delays for different 
excess energies. 

We build on our classical trajectory Monte-Carlo [l^ 
simulation of the classical He*(ls2s) model system, for 
details see 

EEIiM^. We choose He* (ls2s) as a pro- 
totype system to clearly formulate the concept of streak- 
ing two-electron dynamics while avoiding the unneces- 
sary complexity of many-electron systems. However, the 
scheme is not system specific. The same scheme could 
time-resolve, for instance, the collision between Is and 
2s electrons in the ground state of Li |15l-ll7|. Atomic 
units (a.u.) are used throughout this work except where 
otherwise indicated. 

Fig. [1] a) illustrates the concept of the two-electron 
streak camera. First, the Is electron (photo-electron) 
absorbs the XUV attosecond pulse with energies above 
the double ionization threshold. Then, as the electron 
leaves the atom it can collide with the 2s electron and 
transfer some of its energy, resulting in the simultaneous 
ejection of both electrons. The intra-atomic collision is 
typically characterized by the asymptotic inter-electronic 
angle of escape, that can be observed by experiment 
|ll| . To time-resolve the two-electron collision dynamics 
we streak 6^22 t)y adding a weak IR laser field polarized 
along the z axis, F^^(t) = FQ^f(t) cos(cJiRt + 0)z, where (j) 
is the phase between the IR field and the XUV pulse and 
f(t) is the pulse envelope [10]. We choose cjir = 0.0285 au 
(1600 nm) and = 0.007 a.u. (< 1.8 x 10^^ W/cm^) 
so that the streaking field does not alter the attosecond 
collision significantly, but still has an observable effect 
on O12. Here, Fq^ = 0.007 a.u. is chosen to efficiently 
streak excess energies from 10 eV to 60 eV. The IR laser 
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field splits 6>i2(^) in two branches with the lower/upper 
branch corresponding to launching of the photo-electron 
along the ±z direction, see Fig. [T]a). 

Next, we describe how we model the XUV attosecond 
pulse and how its spectral intensity needs to be reflected 
in the weight of trajectories corresponding to different 
excess energies. The electric field of the XUV-pulse is of 
the form: 



p? + pi = 2E- (b) E-«(eV) 



pXUV 



XUV -t'^/4a'^ 



cos(<^^t)i (1) 



with a the standard deviation of the temporal intensity 
envelope I(t). For the current calculations, the spectral 
intensity l{uj) of the XUV-pulse has a FWHM of 75 eV, 
centered at uJq^^ = 120 eV, see Fig.[T]b). The temporal 
intensity envelope I(t) of the transform limited pulse has 
a FWHM of la.u., see inset of Fig. [T]b). In what follows 
we focus on the effect the large energy bandwidth of the 
XUV-pulse has on the streaking process and we neglect 
the effect of the finite FWHM of I(t). The uncertainty of 
the time of photo-absorption will be taken into account 
after the application of the streak camera algorithm as 
an uncertainty in the retrieved delay-times. 

Using first order perturbation theory we compute the 
photo-absorption probability to transition from the ini- 
tial ground state of He*(ls2s) to the final state of double 
electron escape [Tsj : 



Wi^f OC -Cr++(CJ)/(CJ) 



(2) 



with the cross section for double ionization cr++(a;) 
given by (Jabs('^)P^^(<^)- crabs is the cross-section for 
photo-absorption which we calculate in the single elec- 
tron approximation assuming that the Is electron ab- 
sorbs the photon [l9|. P++(co') is the probability for dou- 
ble ionization obtained through our classical simulation 
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[l6| . Finally, we weight each classical trajectory 
for a given photon energy uj by the factor crabs(^)I(^)/^- 
Our goal is to retrieve the delay between photo- 
absorption and ionization of both electrons. Since the 
delay depends on the sharing of the final energy among 
the two electrons [10], we will consider in the following 
only symmetric energy sharing of e < 0.14. The delay 
times for the most symmetric energy sharing correspond 
roughly to the time of minimum approach of the two elec- 
trons, i.e. to the collision time. Here, we have defined the 
energy sharing e = (ei — e2)/(ei + 62), as the dimension- 
less asymmetry parameter between the observable (final) 
kinetic energies 61 and 62 of the two electrons. The analy- 
sis of different energy sharings as described in [III can be 
applied to the following analysis without any restrictions. 
In Fig. [T]c) we plot, a FROG-like trace for two-electron 
ejection, the observable total electron energy as a func- 
tion of (j). In what follows, we describe how we extract 
from Fig. [T] c) the delay times of the intra-atomic two- 
electron collisions for different triggering excess energies. 
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FIG. 1. (a) Sketch of the intra-atomic knock-out mecha- 
nism to be studied with the two-electron streak camera. The 
streaking field causes a decrease in when the photo-electron 
is launched along the H-z direction, since adding ApiR to each 
of the electron momenta results in O12 < ^12- (b) Spectral in- 
tensity of the XUV-pulse scaled by (F^^^)^ Dotted curve: 
Wi^f in the presence of the XUV and IR field averaged over 
all (jy's. (c) Observable E^^^ total electron energy as a function 
of considering excess energies from 4 eV to 60 eV in steps of 
2 eV and double ionization events corresponding to launching 
of the photo-electron (Is) in the zbz direction for (f) ranging 
from 0° to 180°. To illustrate the difference between launch- 
ing of the photo-electron in the + z versus the -z direction 
we plot the E^*^ corresponding to + z for cj) ranging from 0° 
to 180° and the E^*^ corresponding to - z for ranging from 
180° to 360°. 



We first study the effect of the large energy bandwidth 
of the XUV-pulse on streaking the two-electron collision 
dynamics. In Fig. [2] a) we plot the correlation map of 
the excess energy of the XUV-photon and the observ- 
able sum energy E^*^ of both electrons in the presence 
of the streaking IR field. We see that a large range of 
excess energies maps to the same streaked total electron 
energy. Hence, the final electronic state does not cor- 
respond unambiguously to the triggering photon energy. 
For instance, the 20 eV streaked energy maps to excess 
energies ranging from 12 eV to 26 eV. The reason for the 
weak correlation between the streaked and the excess en- 
ergy becomes clear in Fig. [2] c) for 10 eV excess energy: 
the streaked energy changes significantly with <j). 

To retrieve the excess energy from the final electronic 
state with improved accuracy, we introduce a "modified" 
total electron energy, where the effect of the streaking 
IR field is reduced. Therefore, we define a "modified" 
electron momentum pf^°^ by subtracting the momentum 
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FIG. 2. Correlation map of the excess energy E^^ (a) with the 
observable electron energy E^*'^ and (b) with the "modified" 
electron energy E™^'^. The color scale in a) and b) is such that 
the sum of E^^ contributing to a certain E^*'^ is normalized to 
100. (c) Streaked electron energy and (d) "modified" electron 
energy as a function of for E""^ = 10 eV excess energy. 



change ApiR due to the streaking IR field, i.e. 



^mod 



(3) 



where the index i = 1,2 labels the two electrons. The 
change in momentum due to the streaking field (neglect- 
ing the Coulomb potential) is given by 



IR 



OJIR 



sin(A0 + 0). 



(4) 



Here, the shift A(j) = cJiRtdeiay with respect to the max- 
imum of the vector potential of the IR field, Atr ((/>), 
(Fig. [1] a) is due to the delayed electron emission [10|. 
Since we want to retrieve tdeiay we set Acj) = when com- 
puting the "modified" electron momentum p[^°^. Hence, 

ApiR ^ 



■ sin (/). Thus, the "modified" energy E"^°^ is 



given by 

E 



E 

i=l,2 



(5) 



Fig. [2] d) shows that the "modified" electron energy 
varies significantly less with cj) compared to the unmod- 
ified, observable energy E^*^ (Fig. [2] c). Consequently, 
gmod -g strongly correlated with the excess energy, see 
Fig. [2]b). The improved correlation at higher excess en- 
ergies is likely due to the faster collision, i.e. the approx- 
imation A0 ?^ is better at higher excess energies. 

We next explain why at = 0°/180° the "modified" 
electron energy is smaller/larger than the correspond- 
ing excess energy for photo-electrons ejected in the -\-z 
direction (it is the other way around for photo-electrons 
ejected in the -z direction). At = 0° each electron expe- 
riences a force from the IR field in the direction opposite 
to the electron's direction of escape. As a result, each 
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FIG. 3. (a) E" 

excess energy, (b) E™°^ for two electrons as a function of ^ 
for 10 eV excess energy; white solid line shows the average 
of the distribution E™""^ in a) times two. c) E""""^ for two 
electrons as a function of for excess energies between 4 and 
14 eV. The white dashed lines enclose the doubly ionizing 



events with 8 



mod 



180 



8 135 
S* 90 



(a) 



10 eV. 



xlO 



(b) 



noo 
^12 



451 





14 


> 


12 








10 


X 






8 




6 



XlO 

I.. 



45 90 135 180 
(|) (degrees) 



180| 



g135 
S* 90 



45 



x10 

\ 

1^12 



241 
22! 

2o| 

X 

H 18 





(d) 



45 90 135 180 

(|) (degrees) 



> 



^mod ^ 20eV 




45 90 135 180 

(j) (degrees) 



^mod _ 20eV 



45 90 135 180 
(j) (degrees) 



FIG. 4. O12 as a function of for "modified" energies around 
lOeV (a) and 20 eV (c) in the presence of the XUV and IR 
field. A(f) is the shift of the maximum of the vector potential 
of the IR field, corresponding to a maximum of the split of O12 
as a function of 0. (b) Excess energies as a function of that 
contribute to the "modified" energy around 10 eV enclosed by 
the white dashed lines in Fig. [3] c) and similarly (not shown) 
for the "modified" energy centered around 20 eV. 



electron slows down and escapes with an energy smaller 
than the electron's final energy in the absence of the IR 
field. In contrast, at = 180° each electron experiences 
a force from the IR field in the same direction as the 
electron's direction of escape. As a result, each electron 
escapes with an energy larger than the electron's final 
energy in the absence of the IR field. To verify that the 
overall change of the "modified" total electron energy 
with (/) is a one-electron effect, we run our simulation 
in the presence of the XUV plus IR field only for the 
photo-electron (the 2s electron is absent). Since for the 
two electron case we only consider equal energy sharing, 
we compare the two electron case for a certain excess 
energy with the one electron case for half that excess en- 



ergy. Indeed, multiplying by two the distribution of the 
"modified" one-electron energy as a function of (j) for an 
excess energy of, for example, 5 eV (Fig. [3] a) and tak- 
ing the average we find that there is excellent agreement 
with the two electron "modified" electron energy at 10 
eV excess energy as a function of 0, see Fig. [3]b). Note 
that in Fig. [3] and in what follows (Fig. |4] b) and d) and 
Fig. [5j) we focus on double ionization events where the 
photo-electron is ejected in the +z direction. 

The increase of the "modified" electron energy with 
(j) for each excess energy E^^ when the photo-electron is 
ejected in the +z direction, forms the basis for the simple 
algorithm we develop to restrict the residual variation be- 
tween excess energy and "modified" electron energy, i.e. 
the range of "modified" electron energies that pertain to 
a certain excess energy for each value of (j). If our algo- 
rithm is dictated mainly by experimentally accessible ob- 
servables, we compute the collision time corresponding to 
E^® by selecting the doubly ionizing events whose "modi- 
fied" electron energy changes from [E^^ — AE/2,E^^] eV 
for = 0° to [E^" + AE/2, E^^] for = 180°. Choosing 
AE to be roughly 8 eV for all excess energies (method 
1), the thus selected doubly ionizing events have excess 
energies close to the desired E^^ for all (/)'s, and we label 
them by f "^o^. Indeed, the double ionization events with 
gmod _ illustrated with the white dashed lines in 

Fig. [3] c), are the events corresponding to excess ener- 
gies close to E^^ = lOeV for all 0's, see Fig. |4] b) and 
similarly for = 20eV see Fig. Hd). We then deter- 

mine the collision time for f = 10eV/20eV in Fig. H 
a) and c) by extracting A0 from the lower branch of the 
inter-electronic angle of escape O12 W We find 

that A(j) is 4.1° /1. 9° corresponding to a collision time of 
2.5 a.u./1.2 a.u. for = 10eV/20eV, respectively. 

We note that the variation in collision time with excess 
energy suggests that the accuracy of A(j) depends criti- 
cally on the resolution in O12. In order to increase the 
robustness of the retrieval algorithm we determine Acj) 
for a range of bin sizes d0i2 = 4° — 9°. We define the av- 
erage of A(j){d0i2) as the collision phase or collision time 
at a given excess energy or f™^^. 

In Fig. [5] a) we show the collision time for f "^^^ rang- 
ing from 10 eV to 56 eV excess energies. We find that 
the collision time decreases with increasing excess energy, 
as expected. Moreover, we find that the agreement be- 
tween the collision times for a certain £^^^ and the cor- 
responding E^* is best for smaller excess energies. The 
reason is that we compute the delay times correspond- 
ing to a certain S^"^^ using AE ^ 8 eV independent of 
the excess energy. This choice of AE describes best the 
rate of increase of the "modified" electron energy with (j) 
for smaller excess energies. However, AE decreases with 
increasing excess energy. As a result, the agreement is 
worse for higher excess energies. 

Accounting for the fact that AE changes with excess 
energy (method 2) we obtain a much better agreement 
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FIG. 5. □ the collision times for "modified" electron energies 
^mod -|^Q gY gy. ^ ^Y^Q collision times for excess 

energies E^^ ranging from 10 eV to 56 eV. Delay time was 
retrieved by a) method 1 and b) method 2. The error bars 
show the uncertainty in the delay times of 0.6 a.u. since we 
change (j) every one degree. 




E'""^ sharing 

FIG. 6. Correlation map of the excess energy with the "mod- 
ified" energy sharing for the doubly ionizing events with equal 
energy sharing in the absence of the IR field. The color scale 
is such that the sum of E™°*^ energy sharings contributing to 
the equal energy sharing double ionization events for a certain 
E""^ is normalized to 100. 



between the two sets of collision time, see Fig.[5]b). One 
way to do so is by labelling as S^"^^ the doubly ionizing 
trajectories with E"^°^ within ± 1 eV of twice the aver- 
age E™°^ for the one electron problem, see Fig. [3] a) and 
b). However, in this latter algorithm we use as input the 
computed E"^°^ for one electron for each excess energy 
while in the previous paragraph we only use experimen- 
tally accessible data. Finally let us note that both in 
Fig. [5] a) and b) the collision times are computed for 
equal "modified" energy sharing. Similarly to the total 
"modified" electron energy which is strongly correlated 
to the excess energy (Fig. [2] b) the "modified" equal en- 
ergy sharing is strongly correlated to the equal energy 
sharing in the absence of the IR field (Fig. [6]); it is the 
delay times corresponding to the latter energy sharings 
that we aim to streak. 

Concluding, we have demonstrated that the two- 
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electron streak camera can be experimentally realized. 
We obtain a complete picture of the single-photon dou- 
ble ionization process by computing the intra-atomic 
two-electron collision times for different excess energies. 
While our study has been performed in the context of 
two electron escape in atoms it opens the way for time- 
resolving collision dynamics during multi-electron escape 
in atomic and molecular systems. 
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